BUTP-97/28 



Instanton classical solutions 
of SU(3) fixed point actions on open lattices!^ 



Federico Farchioni and Alessandro Papa 
Institute for Theoretical Physics 
University of Bern 
Sidlerstrasse 5, CH-3012 Bern, Switzerland 



February 2008 



Abstract 

We construct instanton-like classical solutions of the fixed point action 
of a suitable renormalization group transformation for the SU(3) lattice 
gauge theory. The problem of the non-existence of one-instantons on a 
lattice with periodic boundary conditions is circumvented by working on 
open lattices. We consider instanton solutions for values of the size (0.6- 
1.9 in lattice units) which are relevant when studying the SU(3) topology 
on coarse lattices using fixed point actions. We show how these instanton 
configurations on open lattices can be taken into account when determin- 
ing a few-couplings parametrization of the fixed point action. 
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1 Introduction 



Fixed point (FP) actions [Q are classically perfect lattice actions, in the sense 
that their spectral properties are free of cut-off effects at the classical level. 
Moreover, they possess scale-invariant instanton solutions up to a minimum 
size of the order of one lattice spacing [0, |^. The latter property, in conjunction 
with a suitable definition of the topological charge operator, allows a theoret- 
ically sound approach to the topology on the lattice |§, ^. In particular, 
since the FP action takes the continuum value on the classical solutions, the 
instanton size distribution in the canonical ensemble is just that of the contin- 
uum up the cut-off scale. In the case of the Wilson action, on the contrary, the 
instanton action is subject to a very fast decrease for decreasing size, and, as 
a consequence, short-ranged topological fluctuations are over-produced in the 
thermal ensemble. These cut-off-scale configurations, which are in fact lattice 
artifacts, hamper the accurate evaluation of the topological quantities (e.g. the 
topological susceptibility x) by Monte Carlo simulations, since their topological 
charge cannot be unambiguously defined on the lattice. 

For a given converging renormalization group (RG) transformation, the FP 
action of any lattice configuration is well-defined and can be determined numer- 
ically with arbitrary accuracy by multigrid minimization Q. In Monte Carlo 
simulations of SU(3) lattice gauge theory, however, due to computing time limi- 
tations only few-couplings parametrizations of the FP action can be used. These 
parametrizations are obtained by a fit procedure on the values of the FP action 
of a representative set of thermally generated lattice configurations Q . Even if 
for a parametrized form the theoretical properties of the FP action hold only in 
approximate sense, parametrized FP actions have shown practically no cut-off 
dependence in Monte Carlo simulations for spectral quantities like the string 
tension and the free energy density even on very coarse lattices ^. 

Topology deserves a separate discussion, since the implementation in a pa- 
rametrized version of the properties of the (exact) FP action concerning topo- 
logical classical solutions introduces additional difficulties. Indeed, classical so- 
lutions never appear among the thermal configurations of the equilibrium en- 
semble, which are the kind of configurations considered in the above mentioned 
parametrizations; if one wants to include the information of the scale-invariance 
of the FP action in the parametrization procedure, classical solutions of the FP 
action must be "artificially" constructed and included in the set of configura- 
tions considered in the fit. In consideration of this, it is not a surprise that 
the first parametrizations adopted in Monte Carlo simulations, e.g. the "type 
Ilia" action proposed in Ref. 0, which perform amazingly well in the case of 
spectral properties, show (as we will see in the following) practically no special 
properties for topology when compared to the standard discretization, i.e. the 
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Wilson action. 

It has been observed ^, ^ that the parametrized FP actions exhibit physical 
scaling already for lattices with spacing a < 1 fm. Since this is just the scale 
of lengths where the instanton size distribution of the continuum reaches its 
maximum the typical instanton size coming into play for such lattices is 
of the order of the lattice spacing. The consequence is that, if the topological 
properties of the continuum theory have to be fulfilled in this precocious scaling 
region, one should concentrate on instantons of size in lattice units /5 ^ 1^ 

The aim of this work is to construct, for this relevant region of size, instanton- 
like classical solutions of the FP action of the type III RG transformation |Q. 
These lattice configurations can be subsequently included in the fit procedure 
for the determination of a new action parametrization. This parametrized FP 
action is expected to well reproduce the FP action of both thermal lattice config- 
urations and classical solutions^ The choice of the type III RG transformation 
is motivated by its property to bring to a FP action with a quite short interac- 
tion range and therefore suitable to be parametrized in terms of a few loops 
extended over 1-2 lattice spacings. 

The main obstacle to this project is represented by the non-existence of 
one-instantons on a lattice with periodic boundary conditions. An effective 
approach in this context, still with periodic b.c, has been proposed in Ref. 
allowing the construction of a few-couplings parametrization of the FP action 
for SU(2) |l^, |ll|. Here, we adopt an alternative procedure, which avoids to 
make use of b.c. at all, i.e. working on open lattices. If on one side this approach 
solves the problem at the root, on the other it introduces many new technical 
subtleties. In this paper we show how all these technicalities can be overcome, 
allowing the construction of instanton classical solutions on open lattices and 
their insertion in the parametrization procedure in a consistent framework. 

An elegant way to accommodate a single instanton on the lattice is also to 
use twisted b.c. Except for the lack of an analytical form for the instanton 
classical solution and for the need to find a proper definition of the instanton 
size, the procedure followed in this paper could be rephrased without problems 
for the case of twisted b.c. 

The remaining of the paper is organized as follows: in Section ^, we briefly 
review the definition of FP action and discuss some consequences of its scale- 

^Here and in the following, the hat on a dimonsionful quantity indicates that it is expressed 
in lattice units. 

''One could deduce that such parametrized FP action in addition should mimic the exact 
FP action even in the case of a topological classical solution with superimposed quantum 
fluctuations. 
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invariance in connection with topology; in Section ^, we present our procedure 
to construct instanton classical solutions of the FP action of the type III RG 
transformation on open lattices; in Section ^, we show how these instanton con- 
figurations can be involved in the parametrization of the FP action; in Section ^, 
we present our conclusions. 



2 Fixed point actions and topology 

The partition function of a SU(N) lattice gauge theory is 

Z = j DU e-P^^''^ , (1) 

where f3A(U) is any lattice regularization of the continuum action expressed in 
terms of products of link variables C/p(ri) G SU(N) along arbitrary closed loops. 
Denoting with /3, C2, C3, . . . the couplings contained in (3A(U), the action can be 
represented as a point in the infinite dimensional space of the couplings. 

A RG transformation in this space of the couplings with scale factor 2 can 
be defined in the following way: 

^-P'A'iV) ^ j j^jj ^-f3[A{U)+T{U,V)] . 

here {V} denotes the coarse configuration living on the lattice with spacing a' = 
2 a, whose links V^(nB) are related to a local average of the original link variables 
U^{n)\ T{U, V) is the blocking kernel which defines this average, normalized in 
order to keep the partition function invariant under the transformation. The 
explicit form of T([/, V) is 

T{U, V)=J2{ Ki^B) - ^Re Tr[ V^inn) QUub) ] ) , (3) 

where Q^i{nB) is a. N x N matrix which represents some average of paths of fine 
link variables U^{n) connecting the sites 2?t,b and 2nB + 2/i; N^{nB) guarantees 
the normalization, while k is a free parameter adjustable for numerical optimiza- 
tion. For the type III RG transformation, (5p(riB) is given by the product of 
the fuzzy links W^(2nB) and W^(271b + /*), each of them being a weighted sum 
of paths of fine links living on the I'' hypercubes which contain also C/^(2nB) 
and Ufj_(2nB + fi), respectively (see Ref. Q for the details, and for the values of 
K and of the other free parameters of the RG kernel). 

On the critical surface /3 = 00, Eq. (J^ leads to the saddle point equation 

A' {V)= mm [A{U)+T{U,V)] ; (4) 
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the FP point of this transformation is therefore 



A'^V) ^mm[A''^{U)+T{U,V)] . (5) 



In the limit /3 — > cx), Af^{nB) is given by 



= ^ max ReTr[M^gt(„B)] • (6) 



A blocking transformation on the link variables Ufj. (n) can be defined, where 
the blocked link Vf_i{nB) is obtained by projecting on SU(3) the average Qfj.{nB)- 
In other words, V^(nB) is given by the SU(3) matrix W maximizing the quantity 

ReTT[WQlinB)] ; (7) 

it is evident from Eqs. ^ and (^) that for a given fine configuration {[/}, 
the corresponding blocked configuration {V} is the solution of the equation 

T{U, V) = at /3 = oo. 

An important consequence of Eq. (|^) is the following: if a configuration {V} 
satisfies the FP classical equations and is a local minimum of A^^{V), then 
the fine configuration {UiV)} minimizing the r.h.s. of Eq. (||) satisfies the FP 
classical equations as well and the value of the action remains unchanged |l^ : 

^FP(T/) = . (8) 

Moreover, since T{U{V),V) — 0, {V} is just the configuration obtained by 
blocking the fine configuration {U{V)}. For this reason, it can be said that 
{U{V)} represents the "inverse-blocked" configuration of {V}. If {V} is an 
instanton-like classical solution of the FP action with size p in coarse lattice 
units, the configuration {U(V)} is also an instanton-like classical solution with 
the same action and size p' — 2p in fine lattice units^ the action is scale- 
invariant. By iteration one can conclude that this scale-invariant instanton 
action is just that of the continuum, ^inst — 47r^/A^. The theorem does not hold 
in the reversed direction: starting from a classical solution, it is not guaranteed 
that the blocked configuration is still a classical solution with halved size in the 
blocked lattice units. In fact, there is a critical instanton size in lattice units 
Pc below which no instanton classical solutions can exist on the lattice even 
with the FP action. This critical size turns out to be typically of the order of 
1 in lattice units Q. In order to check that the blocked configuration of an 
instanton solution is still a solution, it is necessary in principle to verify that 
the inverse-blocking restores the original configuration. 

*The factor two comes from the fact that the size in physical units p does not change under 
blocking. 
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On the basis of the above considerations, a procedure to construct instanton 
classical solutions of the FP action of a given RG transformation up to a size as 
small as the critical one can be outlined as follows on a very fine lattice a 
continuum instanton solution with a large size po is naively discretized; for large 
enough po, the lattice configuration so obtained approximates the corresponding 
classical solution of the FP action with high accuracy. Then, starting from 
this configuration, a certain number of blocking transformations is performed 
in order that the size of the blocked instanton after the last step is of order 
1 in units of the final blocked lattice. As stated in the Introduction, this is 
indeed the relevant region of size when simulating with FP actions. That the 
instanton is not destroyed by the blocking transformation can be checked by 
inverse-blocking, as explained above. 



3 Constructing fixed point classical solutions on 
open lattices 

In this Section we put into practice the procedure outlined at the end of the 
previous Section. The first step is the discretization of a continuum instanton 
solution on the lattice. The size of the instanton in lattice units has to be 
large enough to have a smooth lattice configuration. If one starts from a finite 
lattice with periodic b.c, a difficulty arises since one- instanton solutions with 
such b.c. do not exist in the continuum. One possibility is to take the infinite 
volume one-instanton solution in the continuum, to cut it on a finite lattice and 
to force periodicity. This introduces finite volume effects in the action which 
scale with the lattice size L . This severe volume dependence can be damped 
to order 1 / by performing a singular gauge transformation on the instanton 
configuration at infinite volume, i.e. before cutting it to the finite lattice [p^, [l0|. 
The lattice configuration obtained in this way is the starting configuration of 
the blocking procedure. This is a good approximation of a FP classical solution 
in local sense except in the vicinity of the border, where a systematic deviation 
comes into play. 

Our approach is radically different since it avoids to impose periodic b.c. in 
any step of the full procedure. How this can be done is clear if one considers 
the structure of the blocking kernel T{U,V) (Eq. (^): the blocked link V^(nB) 
is affected only by the fine links entering the definition of Q^i{n^), which, for 
the type III RG transformation, live on the l"* hypercubes of the fine lattice 
containing the fine links U^{2n-Q) and C/p(2nB + fi). Consequently, in order to 
perform a blocking transformation on a lattice {n : = — L/2, . . . ,I//2 — 
1, /i = 1, . . . ,4} without imposing any b.c, it is necessary to know the infinite 
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volume configuration on the larger lattice {n : = — (L + 2)/2, . . . , (L + 2)/2 — 
1, = 1, . . . ,4}. In this way, the blocking transformation is locally equivalent 
to that in the infinite volume. Of course, the blocked configuration (as well as 
the starting one) is known only on a finite sub-lattice of an ideal infinite lattice. 
In the following, we will refer to sub-lattices embedded in an infinite lattice 
without any prescription for their boundary as to "open" lattices^. Given an 
open lattice A with size L and a configuration {[/} living on the infinite lattice 
embedding A, we define {[/'"*} as the set of the links Ufj,{n) of the configuration 
{[/} originating from any site n which belongs to A along both positive and 
negative directions. All the other links of {[/} will be denoted by {If^^^}. 



3.1 Lattice action on open lattices 



Our way of proceeding poses the problem of assigning a lattice action to a config- 
uration living on an open lattice. Even once that a lattice action regularization 
is fixed, there is the further arbitrariness in prescribing which loops are to be 
considered internal and which external to the open lattice. 

We consider lattice action regularizations of the form 

E Q(C)[iV-ReTr(C/c)r , (9) 

C, i>l 

where C denotes any closed path, Uc stands for the product of the link variables 
along the path C and Ci{C) is the coupling associated to the i-th power for the 
loop C (the couplings ci(C) must satisfy the normalization condition which 
ensures the correct continuum limit, see Appendix). For an open lattice A = 
{n : — — L/2, . . . , L/2 — 1, /i = 1, . . . ,4}, after a definite choice of the 
topologies of closed paths C (typically living on a or 2^ hypercube), we 
define the lattice action ^int = X^neA •^{'n) through the clover-averaged action 
density ^ 

An)^^^ E c.(C) [^-^-^-^,^f . (10) 
^ ' N ^ ^ ' perimeter(C) ^ ' 

Here the most important point is that only those loops C which contain the 
site n g A are summed over. The action density A{n) has the right continuum 
limit and when summed on the sites of the infinite lattice reproduces the lattice 
action (|9|). It should be noted that the definition ( p^ ) involves also links which 
are external to the lattice A, but belong to closed paths C containing internal 
sites. Of course, after summing A{n) over the sites of the open lattice A, the 

^ Open lattices in the context of SU (2) topology were used already in Ref . . 
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weights l/perimeter(C) for any given loop topology C sum up to 1 for the 
internal loops, while for the loops on the border they sum to a fractional value. 



If an instanton configuration is known only for the links contained in an 
open lattice, the simplest guess for its full infinite volume action is obtained by 
adding to A-„it the action density of the continuum integrated on the external 
volume^: 



int 



:A- 



■int 



(11) 



where TZ"^ \ A indicates the volume external to the open lattice A. Considering 
the symmetry of the definition of the action density, Eq. (p^), it is natural to take 
as external volume the region identified by a;^ < ~{L + l)a/2 , > (L — l)a/2, 
fj, = 1, ... ,4. Eq. (11) could seem a crude approximation, but as we will show 



shortly it can be quite accurate if some conditions on a/ p and p/L are satisfied. 



In order to check the accuracy of the Eq. 
configurations in SU(2) defined by 



ll), we consider lattice instanton 



[/('"^^t) (n) P exp / (n + sp) ds 



(12) 



where A\j,'^^^\n) are the anti-Hermitian gauge fields for the continuum instanton 
(see Appendix for the explicit expression and for other details needed in the 
following). The definition ([l^ ) corresponds to a lattice instanton obtained by 
an infinite number of blocking steps from the continuum according to the simple 
blocking transformation^ V^i(f^B) — U ^{2n^)U ^{2n^ + p). Expanding in powers 
of the action density ( p^ of the configuration t/^™''*^(n) given in (p^), an 
infinite series is obtained, where the term ~ a^+^" is a linear combination of 
the continuum gauge- invariant operators with nai've dimension 4 + 2?i. The 
leading term of this expansion is just —1/2 Tr( fIiu^*'\x) )^a'*, i.e. the action 
density on the continuum which appears in the integral at the r.h.s. of (pT|). 
For a continuum instanton with size p, this integral depends only on p/L and is 
proportional to {p/L)'^ for small p/L. The first correction to the estimate (nlj) 



is due to the term ~ a in the expansion of the action density (jlCj) and is 
proportional to {a / p)"^ {p / L)^ (see Appendix). Higher dimension operators give 
contributions proportional to (a/p)^"(p/i)^+^", n = 2,3, . . . . The corrections 
are therefore under control when a/p<l and p/L is sufficiently small. In Fig. |] 
we report the estimate (|l^) of the action of a lattice instanton with size p — 8 for 
different values of the size L of the open lattice; we consider three different lattice 
actions, the Wilson action, the plaquette +(1x2) loop tree-level Symanzik 

^Here we adopt the convention of anti-Hermitian gauge fields Af^{x) (see Appendix). 
^Unfortunately, this RG transformation does not converge to a FP in d = 4. 
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Figure 1: Estimate through (|11|) of the infinite volume action of a discretized 
instanton with p = 8 known on open lattices with size L. The action is given 
in units of the continuum value; we consider the Wilson action, the plaquette 
+ (1x2) loop Symanzik action and the plaquette +(2x2) loop Symanzik 
action. The solid lines represent the expected values of the infinite volume lattice 
action: 0.99686923, 0.99998172, 1.00001427, for the three action regularizations 
respectively. 



improved action and the plaquette +(2x2) loop tree-level Symanzik improved 
action^. We observe that the discrepancy between the estimate (11) and the 
theoretical action value goes down rapidly for increasing L. The theoretical 
value is obtained by the expansion of the action in powers of a/ p (see Appendix) 
up to the order 4; for p — 8 this is an excellent approximation. In Fig. |^ we 
perform the same comparison in the less favorable case p = 1.5, for which the 
expansion of action in a/p is not enough accurate. We observe in this case a 
fast convergence to a plateau, corresponding to the exact infinite volume action. 
This value is approximated within some per mill already for L > 6. 



*In the case of tree-level Symanzik improved actions, the corrections to the estimate 
are more suppressed, since they start from (a/p)*(p/L)**. 
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Figure 2: As in Fig. |^ with p = 1.5. In this case a theoretical value for the 
infinite volume action is not available. 



3.2 The blocking procedure 



The next step in the procedure outlined at the end of Section || consists in 
performing repeatedly blocking transformations on smooth instanton configu- 
rations until their size is of the order of 1 in units of the final blocked lattice. 
We denote the instanton configuration on the finest lattice, at the starting point 
of the blocking procedure, by {Uq}; the starting instanton size in lattice units is 
Pq. The instanton configuration after k blocking steps is indicated with {Uk}] 
its size in units of the spacing of the fc-blocked lattice is pk = Po/2'''. 

As we have seen in the previous Subsection, for instantons with size of the 
order of 1 lattice spacing, a 6"* open lattice allows an accurate determination 
of the lattice instanton action at infinite volume through Eq. (pT|). In order 
to calculate on an open lattice the contribution ^int to the total action A 
according to Eq. (pi]), we must be able to calculate also those loops entering 
the definition (10) which lie on the border of the lattice and involve external 
links. This amounts to know the configuration on a lattice larger than the open 
lattice considered. For our purposes, in order to calculate Aint on open lattices 
6^, we need to know the configurations on lattices 10''; this is sufficient, since 
we are interested in loops which do not extend beyond the 2"* hypercube. We 
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have therefore to organize the blocking procedure in order that the lattice size 
after the last blocking step is 10. 

Now, if there were definite b.c, a 10* lattice could be obtained by blocking from 
a 20* lattice. In our approach, where no b.c. are imposed, a 10* lattice can 
be obtained by blocking from 22* lattice, according to what discussed at the 
beginning of Section ||. In the case of three blocking steps, which will turn out to 
be enough for our purposes, the sequence of lattice size is 94 — > 46 — + 22 100. 

The instanton configuration on the finest lattice is obtained by discretization 
according to (|2|) . Since the instanton size halves after each of the three block- 
ing steps, the typical po on the finest configuration has to be around 8 lattice 
units. For such instanton sizes any lattice action gives the same value within 
some per mill (for p = 8, using the a/p expansion up to order {a/ pY, one ob- 
tains 0.99686923 for the Wilson action, 1.00001427 for the tree-level Symanzik 
improved action with the (2 x 2) loop, 0.99998172 for the tree-level Symanzik 
improved action with the (1 x 2) loop). We deduce that such configurations are 
a good approximation of the continuum classical solutions. It must be checked, 
however, that they remain solutions of the FP classical equations of motion with 
the same approximation even after three blocking steps (we recall the discus- 
sion in Section ^) . This can be done by calculating the FP action of the final 
configuration {C/3} and verifying that it reproduces the continuum action. The 
FP action of {C/3} is formally given by Eq. (|^), which in this case reads 

^'^P(C/3)=min[^FP(;7)+r({/,C/3)] . (13) 

One should iterate this equation until the relevant configurations on the finest 
lattice are so smooth that any lattice action can be used in the minimization. 
In practice, for four-dimensional non-Abelian gauge theories, memory and time 
limitations prevent from performing more than one step. Fortunately, after 
one step of inverse-blocking the minimizing configuration {f/min} is already so 
smooth (the action lowers typically by a factor 30-40 Q) that AF^ {U) can be 
well approximated by a tree-level Symanzik improved action; in this context we 
use the plaquette -I- (1 x 2) loop Symanzik action pc| ]. 

Another point to check is that, being {C/3} a FP classical solution, it goes back 
to {C/2} by inverse-blocking. Rigorously, the full procedure should be carried 
out on configurations {C/3} and {C/} living on infinite coarse and fine lattice, 
respectively. Since of course we can work only on finite lattices, it is necessary 
to introduce some simplifications, to be justified within our approach with open 
lattices. 

Our procedure is the following: being {C/3} given on a 6* lattice, we search 
the minimizing configuration {C/min} by locally updating only the links U^(n), 

^We overcome the memory problems associated to such large lattices by allocating smaller 
lattices and exploiting the symmetries of instanton configurations, which we have verified to 
be conserved by the blocking procedure. 
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^j, — 1 , . . . , 4 with n given by 

n = 2nB + ^A<,o- , = 0, 1 , cr=l,...,4 , (14) 

for any site tib of the 6^ coarse lattice. Moreover, we use as trial starting 
configuration in the minimization the configuration {U2} itself. This amounts 
to fix a priori {t/mlnl = {Ui^^} outside the 12^ open lattice whose sites are given 
by Eq. (p^. Since {t/f'^*} is smooth for the values of p/L under consideration, 

it is legitimate to estimate -4ext(C^2) with ^o™"*^ ! moreover, we assume to be 
zero the contribution to T(C/, U^) coming from the part of the summation in 
Eq. (|^) involving sites of the coarse lattice not contained in the open lattice. 

Table 1: Minimized and parametrized FP action in units of the continuum 
instanton action of the {C/3} configurations for several values of p^. The config- 
urations {U^} are known on a 6** open lattice and their infinite volume action 
is estimated as explained in the text. The couplings of the parametrized FP 
action are given in Table Column 2 contains the values of T{Uniin, U3) in 
units of the continuum instanton action resulting from the minimization. 



P3 


T{U„,in,U3)/Anst 






0.1250 


0.007628 


0.0255527 


0.02440536 


0.2500 


0.018860 


0.0552252 


0.05289109 


0.3750 


0.049999 


0.1348518 


0.12942630 


0.5000 


0.100250 


0.2797359 


0.27778473 


0.5625 


0.004224 


0.3737279 


0.38207865 


0.6250 


0.003092 


0.9885723 


0.50021130 


0.6875 


0.002587 


0.9887317 


0.62110525 


0.7500 


0.002215 


0.9893495 


0.73261547 


0.8125 


0.001924 


0.9899657 


0.82569128 


0.8750 


0.001684 


0.9903682 


0.89621836 


1.0000 


0.001306 


0.9903730 


0.97390002 


1.1250 


0.001013 


0.9895449 


0.99372518 


1.2500 


0.000782 


0.9882614 


0.98836493 


1.3750 


0.000602 


0.9867958 


0.97771794 


1.5000 


0.000464 


0.9853111 


0.96955895 


1.6250 


0.000358 


0.9839012 


0.96517944 


1.7500 


0.000278 


0.9826250 


0.96364480 


1.8750 


0.000218 


0.9815158 


0.96376640 



Since the finest configurations {Uq} are built in SU(2) and since the blocking 
is transparent to the color indices, the blocked configurations {U3} are in SU(2) 
as well. On this basis, we searched for {t/min} directly in SU(2), verifying in some 
selected cases that the minimization in SU(3) brings to the same configuration. 
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Coming to the results, we find that for ps larger than a critical size pc — 0.6, 
T{Umin, U3) is of the order of 10^'^, showing that {C/min} is very close to the 
configuration {U2} even in the internal region. Moreover, the quantity 

Ant(f/min)+4xr^+T^(C^mi„,f/3) , (15) 

that is our estimate for reproduces within 1% the continuum ac- 

tion for an instanton (see Fig. g and Table |^). The small deviation is quite 
ps-independent and can be attributed to the fact that the instanton configura- 
tions {Uo} are not exact instanton classical solutions and/or that the tree- level 
Symanzik improved action slightly deviates from the FP action even for the 
smooth configurations of the fine lattice involved in the minimization. 
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Figure 3: Minimized FP action (filled triangles), type Ilia FP action (squares), 
type Illb FP action (diamonds) and Wilson action (circles) of the configurations 
{C/3} in units of the continuum instanton action for different values of the size 
P3. The configurations {U^} are known on a 6* open lattice and their infinite 
volume action is estimated as explained in the text. 



For p3 < Pc the configurations {U3} are not classical solutions anymore since 
A{Uinin)+T{Umin, U3) is niuch bclow the continuum value. These configurations 
{C/3} represent in fact non-topological quantum fluctuations with range one 
lattice spacing. 
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4 Parametrization of the fixed point action 



Once that we have built on a 6* open lattice a set of configurations {U^} rep- 
resenting instanton classical solutions of the type III FP action in the region 
of size yO ~ 1, and have found the corresponding minimizing configurations 
{C^min} — {^^2} in the sense of Eq. (^), we consider in this Section the problem of 
how these configurations can be used to determine a few-couplings parametriza- 
tion of the FP action which keeps track, with a good approximation, of the scale 
properties of the exact FP action. 

Before proceeding in this program, we check how the previous parametriza- 
tions of the type III FP action proposed in Ref. behave on the instanton 
configurations {U^}. These parametrizations, called type Ilia and type Illb, 
were found by fitting the exact FP action, determined by numerical minimiza- 
tion, on a set of ~ 500 thermal Monte Carlo configurations generated on a 2^ 
lattice in the range /3wiison = 5.1 — 50. Both parametrizations are of the form 
and involve the plaquette and the twisted perimeter-6 loop (entries no. 1 and 4 
in Table |[ respectively) with « < 4 Q . In the type Illb parametrization it was 
imposed that the couplings satisfy the on-shell tree-level Symanzik condition at 
the quadratic order in the gauge fields. 

In Fig. ^ we present the behavior of the type Ilia and type Illb actions 
on the configurations {U3} in comparison with the Wilson action. We can 
see that they fail to reproduce the scale-invariant shape of the minimized FP 
action to the same extent as the Wilson action. This is not a surprise since 
these parametrizations were not required to well reproduce the FP action on 
classical solutions, but on thermal configurations. We see also that the tree-level 
Symanzik improvement at the quadratic level in the gauge fields has no effect 
on the action of instanton classical solutions. For the values of /3 where the FP 
actions are expected to display physical scaling 1^, |^ , corresponding to a < 1 
fm, the relevant region in the instanton size is /i ^ 1 in these conditions, 
the behavior exhibited in Fig. || by the type Ilia and type Illb actions could 
be responsible for distortions of the lattice signal when computing topological 
quantities |ll|. Indeed, the small-size (p ~ 1) topological configurations, having 
an action lower than the continuum instanton action, are over-produced in the 
thermal ensemble and their contribution to the total topological signal is over- 
estimated, when a non-zero topological charge is assigned to them |2[ |j . 

Parametrizations of the FP action which approximately possess the scale- 
invariance properties of the exact FP action can be built by taking into account 
the additional information coming from the classical solutions ||]. In the fol- 
lowing we explain in some detail how this has been done in the present case. 
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The fit procedure starts by fixing the form of the parametrization of the FP 
action, Eq. (||) in our case, and, accordingly, a finite number of loop topologies 
and powers. We have restricted ourselves to all the loops living on the 2^ 
hypercube (see Table ||) and to four powers of the trace. 

The first condition of the fit is that the FP equation (^) holds also at a para- 
metrized level 

All.{V)^All.{U^^^)+nu^;^,V) . (16) 

In the framework of the fit procedure this means the minimization of the quan- 
tity 



Ql(Ca) 



E 

{V} 



(17) 



where Cq stands for Ci{C) and T^iU) for 1/N [N — Re Tr([/c')]*, being a a 
collective index for both the loop topology index C and the power index i in 
Eq. (§. 

The second condition is that the exact values of the FP action on the minimizing 
configurations {J/min} are reproduced by the parametrized FP action 



(18) 



the r.h.s. is known only in approximate sense (for the instanton configurations, 
we used the plaquette -I- (1x2) loop tree-level Symanzik improved action). The 
quantity to be minimized is in this case 



(Ca) 



E 

{V} 



Y.cMUn.in{V)) - ^^P(C/,ni„(V^)) 



(19) 



The summation in Eqs. (QTI) and (19) is intended to run on a representa- 



tive enough set of configurations. We have considered a set of 500 thermal 
equilibrium configurationtj^ and, in addition, the instanton configurations {t/3} 
built in the previous Section, retaining the freedom to give different weights to 
the first and to the second conditions of the fit and to thermal and instanton 
configurations. In formulae, the quantity to be minimized has been taken of the 
form 

^ — j^(thorm) I ^(thorm)Q(therm) _|_ j^(therm)Q(thorm) j 

+ A('"^')[ Af"'*^gf"*^ +A^""*^Q^'"''M ■ (20) 

The Q's associated to each sector of the fit are polynomials quadratic in the 
coefhcients Cq,; the quadratic term of each polynomial is characterized by a 

^"The same used in Rof. [|7j; we thank the authors of Ref. |Q for having supplied these 
configurations to us. 
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matrix of order equal to the number of the coefficients Ca entering the action 
parametrization (in our case 112, coming from 28 loop topologies times 4 pow- 
ers). The eigenvalue analysis of these matrices shows that they are characterized 
by all but one quasi-zero modes. For example, in the case of the matrix coming 
out in q'^^^'^™\ the eigenvalues, normalized to give trace 1, are distributed as 
follows: 0.9849 the largest eigenvalue, 0.0138 the second, 0.0007 the third, and 
so on in decreasing order. The presence of so many quasi-zero modes indicates 
that there are many equivalent sets of couplings Cq, satisfying with almost the 
same accuracy a given sector of the fit. In view of this relatively wide free- 
dom, it is reasonable to expect that there exist combinations of the couplings 
Ca which nicely fit all the requested conditions simultaneously. Among different 
equivalent parametrizations, of course the fastest for Monte Carlo simulations 
should be chosen. 



We have proceeded in the following way: first of all, we have restricted 
ourselves to all the possible combinations of no more than three topologies of 
loops among those living on the 2'^ hypercube (see Table ||). We have consid- 
ered combinations containing always the plaquette and one or two other loops; 
moreover, we have excluded parametrizations involving two perimeter-8 loops, 
which are in general too slow for standard Monte Carlo updating algorithms. 
For any chosen combination, we have first performed the fit only on the thermal 
equilibrium configurations, that is we have imposed A^*^''™) = 1 and A^'""*) = 
in Qtot- The weights A^*^*"^™^ and \^^^^™^ have been fixed in order that, for the 
set of couplings minimizing Qtotj both q'^^^^™^ and Q^^^'^™^ are as close as pos- 
sible to their absolute minima, which are realized for A2*^'^"°^) = (1, 0) 
and (^A^*'^''™-', Aj*'^"'^™'*) = (0,1), respectively. We have found many different 
parametrizations giving equally good fits, as expected because of the many 
quasi-zero modes involved in the procedure. We could reproduce in particular 
the parametrization Ilia of Rcf. [Q, which is indeed the most economical for 
numerical simulations. 

In order to fix A^'"'''-' and A2"'^*'' for the same choice of loop topologies, we 
have repeated the procedure considering only the instanton configurations, that 
is for A^*'^°''™^ = and A'^'"*'*-' = 1 in Qtot- Then we have combined thermal 
sector and instanton sector of the fit, starting with A^**^"™) = 1 and A^'"*^*-* = 
and progressively increasing the A^'"''*-'. We have realized that for not too large 
values of A*-™''*-', the quality of the fit in the sector of the thermal configura- 
tions remains good, still a consequence of the presence of many zero modes in 
this sector of the fit. To choose among the different, equally acceptable, loop 
combinations, we have measured each candidate parametrized action on the in- 
stanton configurations {J/s} and checked to what extent it reproduces the values 
of the FP action obtained by minimization (Table ^ . The last test has turned 
out to be the most stringent, leading to a strong restriction of the candidate 
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parametrizations to a few ones involving always the plaquette, a perimeter-6 
loop (entries 2-4 in Table ^ and a perimeter-8 loop (entries 5-28 in Table ||). 
Once a fitting parametrization -4p£. is found, it must be checked that it is pos- 
itive (in the sense that Ap^^ > for any configuration), since this property is 
not a priori guaranteed. We could not find any sufficient condition for pos- 
itivity which was not too restrictive and we had to turn to numerical checks 
on some selected classes of configurations. Of course these checks, although 
quite stringent, do not prove in absolute the positivity of the action. We found 
that the positivity test represented a real bottleneck in the selection of the 
parametrizations, since as soon as instantons were introduced into the fit the 
candidate parametrizations lost positivity in almost all the cases. We could 
escape this problem only in one case, namely for the parametrization given in 
Table || which involves the plaquette, the bent rectangle (entry no. 3 in Table ||) 
and the twisted perimeter-8 loop (entry no. 25 in Table ||). Of course, the simu- 
lation time cost of this parametrization is high, since it involves loops with large 
multiplicity and large perimeter. This action is 5 times slower than the type 
Ilia parametrization of Ref. and 35 times slower than the Wilson action with 
the same Monte Carlo algorithm. In Fig. ^ and in Table [| this parametrized 
action is compared with the exact FP action on the configurations {C/3}. It can 
be seen that the agreement is good for instanton size ps > 0.8. In the region 
0.5 < ps < 0.8 the quality of the fit is poor: actually we excluded from the 
beginning the instanton configurations corresponding to this region, since here 
the Symanzik tree-level improved action used in the minimization procedure 
to determine the FP action is not completely reliable for the smallness of the 
involved sizes. For < 0.5 the agreement is again very good: since below the 
critical size the blocked configurations represent small range fluctuations around 
the trivial vacuum, they resemble thermal configurations; so, the agreement in 
this region is actually an indication of the quality of the fit in the sector of the 
thermal configurations. 



A few-couplings parametrization of the FP action well reproducing the in- 
stanton action was already found in |11 using periodic b.c. according to the 
technique described in Also in this case it was found that an additional 
perimeter-8 loop is needed to well parametrize the FP action in presence of 
instanton solutions, although a different choice of closed paths turned out to be 
preferable, i.e. the loops no. 1, 4 and 19 of Table ||; the computation time of 
this parametrization is roughly the double of that of Table ||, the behavior of 
the action with the instanton size being qualitatively the same. 



In order to determine topological quantities from MC simulations a further 
step is needed, namely a suitable definition of the topological charge operator. 
The use of the FP topological charge operator Qpp in conjunction with the FP 
action allows a consistent definition of the topology on the lattice ^, ||, |[ |l^, |ll| , 
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since in this case all the configurations with non-zero topological charge have ac- 
tion larger or equal than the continuum instanton action in the given topological 
sector [Q. In particular, this implies that there is no over-production of small- 
size instanton configurations. This ideal picture is not automatically realized 
when a parametrized version of a FP action is used in numerical simulations. In 
this case it could well be that there exist configurations having (parametrized) 
FP action lower than the continuum instanton action and topological charge 
different from zero. Our results (Fig. U) seem to suggest that the FP topo- 
logical charge and action reproduce the continuum values up to size p ~ 0.5, 
while the parametrization of Table || deviates from the ideal behavior at larger 
values (around 0.8). This indicates that further improvements are needed: one 
possibility is to improve the parametrization itself by adding new loops or by 
investigating other functional forms of the parametrized action; another is to 
improve the inverse blocking procedure by increasing the number of steps or by 
using a better approximation of the FP action in the minimization procedure. 
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Figure 4: As in Fig. ^ with the minimized FP action (filled triangles) and the 
parametrized FP action defined in Table ^ (squares). 



We turn now to some questions concerning the Symanzik improvement in 
the context of the parametrization of a FP action. For the technical details 
involved, we refer to the Appendix. 

The FP action satisfies by definition the on-shell tree-level Symanzik condi- 
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tions to all order in a} . For instanton configurations the Symanzik conditions at 
order a" ensure the absence of the term 0{a/ p)^ in the expansion of the action 
in powers of a/ p. Of course, a parametrized form of the FP action does not nec- 
essarily fulfill the Symanzik conditions. Nevertheless, a parametrized form of 
the FP action having exactly vanishing cut-off effects at tree-level at the lowest 
order in is preferable, since it would be automatically improved in the per- 
turbative regime and everywhere smooth configurations are physically relevant. 
The on-shell tree-level Symanzik conditions at order a? would in particular en- 
sure that the first correction to the continuum action of instanton solutions is 
0{a/ pY; as a consequence, the behavior of the parametrized action on instan- 
tons with large size (not included in the fit procedure) would be improved. In 
consideration of this, we have repeated the fit procedure including with some 
weight the on-shell tree- level Symanzik conditions at order a? . However we were 
not able to find any positive definite parametrization exhibiting a nice fit and 
supporting also the Symanzik conditions. 

We conclude observing that the parametrization of Table ^ badly violates 
the Symanzik conditions. This is not in contradiction with the fact that this 
parametrization nicely fits the FP action for the restricted set of configurations 
taken into account in the fit procedure, since a compensation effect between 
contributions of different orders in (equally important for the configurations 
considered) occurs. 



5 Conclusions 



The conclusion of this work is that it is possible to carry out in a consistent way 
the program of building instanton classical solutions of the FP action of a RG 
transformation for the SU(3) lattice gauge theory without imposing any b.c. We 
have shown how these configurations, defined on open lattices, can be included 
in the fit procedure to determine a few-couplings parametrization of the FP 
action which reproduces the scale-invariance properties of the exact FP action 
up to a size p as small as 0.8. Some additional work is needed to improve the 
parametrization in order to make it effective in numerical MC determinations 
of topological quantities. 
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Table 2: List of the loops living on a 2^ hypercube: they are classified according 
to the 28 different topologies. For each kind of loop, a representative orientation, 
the multiplicity and the contribution to the operators Oq, Oi, O2 and O3 (see 
Appendix) are given. 



loop 


path 


multiplicity 




ri 


r2 


rs 


1 


X, y, -X, -y 


6 


1 


1 








2 


X, y, y, -X, -y, -y 


18 


8 


20 








3 


X, y, z, -y, -X, -z 


72 


16 


4 





12 


4 


X, y, z, -X, -y, -z 


24 


8 


-4 


4 


4 


5 


X, y, X, y, -X, -y, -X, -y 


24 


8 


32 








6 


X, y, X, y, -X, -X, -y, -y 


48 


36 


132 








7 


X, y, X, z, -y, -X, -X, -z 


384 


192 


48 


96 


144 


8 


X, y, X, z, -X, -z, -X, -y 


192 


32 


8 





24 


9 


X, y, X, z, -X, -y, -X, -z 


96 


40 


16 


24 


24 


10 


X, y, X, z, -X, -X, -z, -y 


384 


160 


256 





96 


11 


X, y, X, -y, -X, y, -X, -y 


24 





-12 








12 


X, y, y, X, -y, -X, -X, -y 


24 





-24 








13 


X, y, y, z, -y, -z, -X, -y 


192 


32 


56 





-24 


14 


X, y, y, z, -y, -y, -X, -z 


192 


80 


80 





96 


15 


X, y, y, z, -y, -X, -z, -y 


192 


48 


24 


-24 


72 


16 


X, y, y, z, -X, -y, -y, -z 


96 


72 


24 


48 


48 


17 


X, y, y, -X, z, -y, -y, -z 


96 


64 


112 





48 


18 


X, y, y, -X, -X, -y, -y, X 


12 


16 


64 








19 


X, y, z, t, -z, -t, -X, -y 


192 


32 


32 








20 


X, y, z, t, -z, -y, -X, -t 


192 


48 


-24 





72 


21 


X, y, z, t, -z, -X, -t, -y 


384 


96 








96 


22 


X, y, z, t, -z, -X, -y, -t 


384 


128 


-112 


48 


144 


23 


X, y, z, t, -y, -X, -t, -z 


96 


32 


-16 





48 


24 


X, y, z, t, -y, -X, -z, -t 


192 


80 


-112 


48 


96 


25 


X, y, z, t, -X, -y, -z, -t 


48 


24 


-48 


24 


24 


26 


X, y, z, -y, -z, y, -X, -y 


96 


16 


28 





-12 


27 


X, y, z, -y, -X, y, -z, -y 


48 


4 


16 


-12 


12 


28 


X, y, -X, -y, X, y, -X, -y 


12 


4 


4 









Table 3: Couplings of the parametrization of the type III FP action determined 
in this work (see definition (H)). 



loop 


Cl 


C2 


C3 


C4 


plaquette 


-0.544089 


1.940530 


-0.466298 


0.016159 


bent rectangle 


0.109877 


-0.114994 


-0.008653 


0.007343 


twisted perimeter-8 


-0.008914 


0.004823 


0.015426 


-0.001432 
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A Appendix 



We consider lattice instanton configurations in SU(2) defined by ||T^ 

U^{x)=Pexpf Af,{x + sil)ds , (A.l) 
Jo 

where Afj_{x) = 'iJ2a^'^^li^)'^'^/'^ ^^'^ anti-Hermitian gauge fields of the infi- 
nite volume continuum instanton solution (cr°, a — 1,2,3, are the Pauli matri- 
ces). They are given by 



(A.2) 



where rj'^^, a = 1,2,3, is the 't Hooft tensor p is the instanton size and 
we have assumed that the instanton is centered at the origin. The field tensor 
Ff_,i, = d^A^ - d^Af^ + [A^j^,Atj] turns out to be 



-t- p^y 



(A.3) 



A lattice action A — A{x) defined in terms of the symmetrized action 
density 

, 1 , „, [Af - Re Tr(C/c)]^ 

Ax) = ^ 2^ c,(c)- 



C3x,i>l 



pcrimeter(C) 



(A.4) 



can be expanded by using (^) as ^ = Y.'^^^ a^+'^'^ 0^^+^'^\ where 0(4+2") 
is a combination of all the independent gauge-invariant operators with nai've 
dimension A + 2n. The coefficients of this combination depend on the coefficients 
Ci{C) for a given choice of the topology of the loops C involved in the definition 
of the action. The expansion of the action up to the order is 



Y,r^{C)c^{C)Y,0^{x) 



L C 



MC)c,{C) J2 02{x) + ^ r3(C)ci(C) 0,{x) 



with 



Ooix) 
Oi{x) 



(A.5) 

(A.6) 
(A.7) 
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02{X) 



(A.; 



(A.9) 



being V^F^x = d^F^\ + [A^, Fj^a]- The coefReients rg, ri, r2 and are given in 
Table ||for all the loops whi ch li ve on a hypercube 2'*. They have been found by 
considering the expansion (A. 5) on configurations representing constant fields 
with low intensity. The normalization condition ''o(C')ci(C) — 1 has to be 
satisfied in order to have the correct continuum limit. The on-shell tree-level 
Symanzik improvement at O(a^) is realized by imposing X^c ~ ^ 
and Xp 7'2(C)ci(C) = 0; no condition comes from the operator 03(0;), since it 
is identically zero on configurations which satisfy the equations of motion 118]. 



On an instanton classical solution (A. 2), the operator O'iix) is identically 
zero, while the other operators are 



Co,inst(2;) 



(a;2 + p2)4 
-48p4a;2 

-192p4a;2 
(a;2 + p2)6 



(A.IO) 



We observe that the operators Oi{x) and 02(2;) are not independent when 
acting on the instanton configurations. On such configurations, the O(a^) term 
in the action vanishes if the weaker condition Y^dri^C) + 4r2(C))ci(C) = is 
satisfied. The contribution of the operators 01(0;) and 02{x) to the action ( |A.5|) 
inside a hypersphere with radius R centered in the origin can be estimated by 
integrating the expressions on the r.h.s. of Eqs. ( A.IO ) on the continuum. It 
turns out that 



d'^x Oo,i 



|2;|<i?. 



d X Olinst 



|a:|<-R 



d X 02.^81 



87r2 



87r2 



Sir' 



87r2 



4 a' 



1 



l + 3(i?/p)^ 

-1]3 

O 

1 - 



'|2:|<_R 



m/py 
p_ 

R 
a 
P 
a 
P 

d'^x Oi 



[{r/pY + If 



(A.ll) 



|2;|<_R 
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It can be seen from these results that the integration of an operator with nai've 
dimension 4 + 2n outside a region with characteristic size R in the Euchdean 
four-dimensional space scales as (a/p)^"(/9/i?)^+^". 



Taking the infinite volume limit in Eqs. (A.ll), the lattice action of an 
instanton with size p can be estimated as 



A 



8^2 



^ (^) ' E( ^i(^) + ) ci(C) + O Q ' 



(A.12) 



For the Wilson action, the plaquette + (1x2) loop tree-level Symanzik improved 
action (ci(plaq) ~ 5/3, ci(l x 2) = —1/12) and the plaquette +(2x2) loop 
tree-level Symanzik improved action (ci (p laq) = 4/3, ci(2 x 2) = —1/48), also 
the terms 0(a/p)* of the expansion ( A.12| ) are known |l6| : 



A 



Wilson 



1-; - - 



A: 



Sym,lx2 



70 \p 

17 
210 \p 



As 



ym,2x2 



^"^ 35 U 



01^ 



O 



(A.13) 
(A.14) 
(A.15) 
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